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About 1/3 of X-ray-luminous clusters show smooth, unpolarized radio emission on ~Mpc scales, 
known as giant radio halos. One promising model for radio halos is Fermi II acceleration of seed 
relativistic electrons by turbulence of the intracluster medium (ICM); Coulomb losses prohibit accel¬ 
eration from the thermal pool. However, the origin of seed electrons has never been fully explored. 

Here, we integrate the Fokker-Planck equation of the cosmic ray (CR) electron and proton distri¬ 
butions in a cosmological simulations of cluster formation. For standard assumptions, structure 
formation shocks lead to a seed electron population which produces too centrally concentrated radio 
emission. Instead, we present three realistic scenarios that each can reproduce the spatially flat radio 
emission observed in the Coma cluster: (1) the ratio of injected turbulent energy density to thermal 
energy density increase significantly with radius, as seen in cosmological simulations. This generates 
a flat radio profile even if the seed population of CRs is steep with radius. (2) Self-confinement of 
energetic CR protons can be inefficient, and CR protons may stream at the Alfven speed to the 
cluster outskirts when the ICM is relatively quiescent. A spatially flat CR proton distribution de¬ 
velops and produces the required population of secondary seed electrons. (3) The CR proton to 
electron acceleration efficiency K ep ~ 0.1 is assumed to be larger than in our Galaxy ( K ep ~ 10 —2 ), 
due to the magnetic geometry at the shock. The resulting primary electron population dominates. 

Due to their weaker density dependence compared to secondary electrons, these primaries can also 
reproduce radio observations. These competing non-trivial solutions provide incisive probes of non 
thermal processes in the high-/3 ICM. 

PACS numbers: 

About one third of X-ray-luminous clusters show 
smooth, unpolarized radio emission on ~Mpc scales, 
known as giant radio halos (RHs) [TJ. They appear only 
in disturbed, merging clusters and the RH luminosity 
correlates with the X-ray luminosity [2,3 and the Comp¬ 
ton y-parameter jij. The RHs show that CRs and mag¬ 
netic fields permeate a large volume fraction of the intra¬ 
cluster medium (ICM). The dominant CR source, given 
the smoothness and enormous extent of RHs, is thought 
to be structure formation shocks @,0|- At the same time, 
plasma processes, the origin of magnetic fields and par¬ 
ticle acceleration in a turbulent, high-/3 plasma like the 
ICM are not well understood. Radio halos thus provide 
an incisive probe of non-thermal processes in the ICM. 

One promising model for RHs is re-energetization of 
seed suprathermal electrons by Fermi II acceleration 
when ICM turbulence becomes transonic during merg¬ 
ers Due to the short radiative cooling time of 

high-energy relativistic electrons, the cluster synchrotron 
emission quickly fades away after a merger, which natu¬ 
rally explains the observed bimodality of RHs. 

However, there is a salient piece missing in the tur¬ 
bulent reacceleration model. It relies heavily on the as¬ 
sumption of an abundant, volume-filling population of 


seed suprathermal electrons; direct Fermi II acceleration 
from the thermal pool is precluded by strong Coulomb 
losses 30- These seeds are presumed to be either fos¬ 
sil CR electrons (CRes) accelerated by diffusive shock 
acceleration (DSA) during structure formation 11], or 
secondaries injected by hadronic interaction of CR pro¬ 
tons (CRps) with thermal protons [12] ■ While analytic 
estimates have been made, there has been no ab initio 
demonstration that structure formation can lead to the 
required abundance of seed electrons with the correct 
spatial and spectral characteristics. This is a non-trivial 
requirement: Coulomb cooling in dense cluster cores is 
severe, and DSA fossil electrons may not survive. On the 
other hand, for secondaries to constitute the seed popu¬ 
lation, the CRp population required in the best-studied 
case of the Coma cluster must have a very broad and flat 
(or even slightly inverted) spatial profile [12] , in contrast 
with the thermal plasma whose energy density declines 
steeply with radius. In this Letter we show that such a 
distribution is not predicted by cosmological simulations 
(see lower right panel of Fig. [T]) [see also 0E3. 


Indeed, arriving at a seed population with the required 
characteristics is highly constraining, and has the poten¬ 
tial to teach us much about the origin of CRps/CRes in 
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clusters. We consider 3 new possibilities: (i) Our model 
M-turbulence: a significantly flatter turbulent profile 
than what was adopted in |13[, which allows seed CRps 
to follow the steep profile that is suggested by struc¬ 
ture formation simulations, (ii) Our model M-streaming: 
streaming CRps that produce flat distributions of CRs 
in the ICM 0, 0, which also flattens the secondary 
electron distribution, (iii) Our model M-primaries: if 
the acceleration efficiency of CRps is below about 0.1 % 
in weak (perpendicular) shocks and the ratio of injected 
electrons-to-protons K ep ~ 0.1, this yields a dominant 
primary population with a flat spatial distribution, since 
primaries have a weaker density dependence than secon¬ 
daries. In this work we pursue these three possibilities 
further. We employ cosmological simulations of CRs in 
clusters, in tandem with new insights from our recent 
work on DSA generated fossil electrons [18J, to gener¬ 
ate the first quantitative calculation of primary and sec¬ 
ondary seed electrons. 

Method. The transport of relativistic electrons and 
protons in the ICM is a complex process that depends 
both on the details of the thermal component (gas den¬ 
sity, temperature, and pressure) as well as non-thermal 
component (turbulence, magnetic fields, fossil CRs). We 
use high resolution galaxy cluster simulations to derive 
the thermal and fossil CR properties (shock accelerated 
primary CRes and CRps, as well as secondary CRes pro¬ 
duced in p-p collisions) p, a! [3, 0 ■ hr this Letter we 
focus on our simulated cluster, g72a, which is a massive 
1.6 x 10 15 -Mq cluster that experienced a merger about 
1-2 Gyrs ago. Since the cluster mass, density and tem- 
peratureprofiles are all similar to the well studied Coma 
cluster 00 , we will compare our calculations to radio 
and gamma-ray observations of Coma. 

In our Galaxy, the CRe-to-CRp ratio at a few GeV is 
K ep ~ 10 -2 . Hence, we adopt this as a fiducial value for 
the CRe-to-CRp acceleration efficiency (see 0 for more 
discussion). However, as recent PIC simulations have 
shown, this is likely very different at weak shocks, with 
electrons efficiently accelerated at perpendicular shocks 
211 ] , and ions efficiently accelerated at parallel shocks 



221 ]. Thus, depending on magnetic geometry, K ep could 
be either larger or smaller. In this work we use a sim¬ 
ple test-particle model for the CRp acceleration 18, 23j . 
The ratio of accelerated proton-to-dissipated energy in 
the downstream of strong shocks varies from 1 - 10 %, de¬ 
pending on the adopted model (for more details, see the 
Results section), and is a factor 10-100 lower for weak 
shocks. However, some observations of radio relics sug¬ 
gest higher values of K ep , due to the absence of gamma- 
ray emission, which probes the CRp population [24]. 
This suggests primary CRes as a viable alternative sce¬ 
nario to secondary CRes as seeds for the giant RHs. In 
our M-primaries scenario, we adopt K ep = 0.1 (viable for 
primarily perpendicular shocks) to test this possibility. 

As previously noted, secondaries produced by shock ac¬ 


celerated CRp have the wrong spatial profile to explain 
RH observations; because they arise from a two body 
process, they are too centrally concentrated. They also 
produce 7 -ray emission in excess of Fermi-LAT upper 
limits 000 . However, if CRps stream in the ICM, 


then their spatial profile could potentially flatten suffi¬ 
ciently 0 17 1 . This scenario is very attractive: it gen¬ 
erates seed electrons with the right spatial footprint, and 
by removing CRps from the core, obeys gamma-ray con¬ 
straints. Turbulence plays two opposing roles: Alfvenic 
turbulence damps waves generated by the CR stream¬ 
ing instability 0 , 0 ], thus reducing self-confinement; 
but compressible fast modes scatter CRs directly. Tur¬ 
bulent damping is still efficient for highly subsonic con¬ 


ditions [17]. while we assume compressible fast modes 


only provide effective spatial confinement during the pe¬ 
riods of transonic, highly super-Alfvenic (Ma ~ 5) tur¬ 
bulence associated with mergers. Thus, CRs can stream 
out when the cluster is kinematically quiescent. Further¬ 
more, even Alfvenic streaming timescales are relatively 
short (~ 0.1 - 0.5 Gyr; 0 ) compared to the timescale 
on which the CRp population is built up. Based on these 
findings, we adopt a toy model for our M-streaming sce¬ 
nario in which CR streaming quickly produces flat CRp 
profiles. We assume that CRs cannot stream significantly 
past perpendicular R-fields at the accretion shock, so 
that the total number of CRs is conserved. 

Given a seed population of CRs, we adopt essen¬ 
tially the same set of plasma physics assumptions as 
the reacceleration model for RHs [1, 0 , leaving explo¬ 
ration of parameter space to future work. We solve the 
isotropic, gyro-phase averaged Fokker-Planck equation 
(via a Crank-Nicholson scheme) for the time evolution 
of the CRe distribution in the Lagrangian frame 00 : 
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Here / e is the one-dimensional distribution in position x 
(suppressed for clarity), momentum p and time t (which 
is normalized such that the number density is given by 
n e (t) = f dpf e (p , t)), d/dt = d/dt+v-V is the Lagrangian 
derivative, v is the gas velocity, \dp/dt\ represents radia¬ 
tive (r) and Coulomb (C) losses, D pp is the momentum 
space diffusion coefficient, and Q e denotes the injection 
rate of primary and secondary electrons in the ICM. The 
V • v terms represent adiabatic gains and losses. During 
post-processing of our Coma-like cluster simulation, we 
solve the Fokker-Planck equation over a redshift inter¬ 
val from z = 5 to 0. The simulated cluster undergoes 
a major merger over the last 1-2 Gyrs that injects large 
turbulent eddies. After about 1 Gyr those have decayed 
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down to the scale needed to reaccelerate particles. In 
all our calculations we assume that turbulent reaccel¬ 
eration is efficiently accelerating particles for 650 Myrs 
and that during this turbulent phase CR streaming and 
spatial diffusion can be neglected. In our M-streaming 
model, CR streaming and diffusion are incorporated sep¬ 
arately during kinematically quiescent times that precede 
the merger. As a result, flat CRp profiles are produced 
on relatively short timescales (~ 0.1 — 0.5 Gyr). 

The time evolution of the spectral energy distribution 
of CRps, / p , is similarly given by: 


dp 

dt 
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where Q p (p, t) denotes the injection rate of shock acceler¬ 
ated CRps as a function of momentum p and time t, and 
Thad is the timescale of hadronic losses that produce pions 
via CRp collisions with thermal protons of the ICM [e.g. 
[l2| . We incorporate momentum diffusion for electrons 
and protons from transit-time-damping (TTD) resonance 
with compressible magneto-hydrodynamic (MHD) tur¬ 
bulence, to model Fermi-II reacceleration || T2|. The 
TTD resonance requires the wave frequency to = knv\\, 
where k\\ and U|| are the parallel (projected along the 
magnetic field) wavenumber and particle velocity, re¬ 
spectively. This implies that the particle transit time 
across the confining wave region matches the wave pe¬ 
riod, A||/u|| = T. The resonance changes the component 
of particle momentum parallel to seed magnetic fields, 
which over time leads to increasing anisotropy in the par¬ 
ticle distribution that decreases the efficiency of reaccel¬ 
eration with time. As in ref. fl 2 | , we assume that there 
exists a mechanism—such as the firehose instability— 
that isotropizes the CR distribution function at the gy- 
roscale and on the reacceleration time scale, which en¬ 
sures sustained efficient reacceleration with time. The 
particle pitch-angle averaged momentum-diffusion coef¬ 
ficient of isotropic particles that couple to fast magne- 
tosonic modes via TTD resonance is [ 8 j (Eqn. 47): 


1 9 = 


rarccos(V ph /c) s j n 3 g 

J 0 |cos 9 1 



Here Vj,h is 


the phase velocity of the fast magnetosonic waves given 
approximately by the sound speed, Eph ~ c s . For a sound 
speed typical for the ICM of 1000 km/s, Ig ~ 5. As 
in [§J, we initially assume that the velocity of turbulent 
eddies is Vq ks 0.47c s throughout the cluster. This gives 
a turbulent acceleration time scale, t pp = p 2 /AD PP , that 
is typically few 100 Myrs in the ICM. 

We adopt a simplified isotropic Kraichnan MHD tur¬ 
bulent spectrum for the fast modes per elemental range 
dk of the form 


W(fc)« ^I 0 p(Vp h )k~ 3 /\ 


(4) 


for k 0 < k < k cu t, where we assume an injection scale 
for the turbulence, fco = 27r/(100 kpc). The volumet¬ 
ric injection rate of turbulent energy, Jo, is fixed by re¬ 
quiring that the total turbulent energy density on the 
largest scales e t urb = fW(k)dk ~ 0 . 2 e t h, where e t h is 
the thermal energy density |sj, [l2] • In this work we in¬ 
vestigate different spatial models for injected turbulence. 
We assume that e tU rb ex e“j) u , where a tu = 0.69 for M- 
turbulence , atu = 0.84 for M-streaming , and «t u = 0.91 
for M-primaries (note that in previous work, cnt u = 1 
was adopted E3)- Our flatter turbulent profiles are mo¬ 
tivated by fits to cosmological simulations 29j-31| and 
the range indicates uncertainties of the turbulent profile 
in Coma. Future observations (by Astro-H) and sim¬ 
ulations will help to clarify this issue. Provided dissi¬ 
pation of turbulence in the ICM is collisionless, turbu¬ 
lent cascades of compressible modes become suppressed 
when thermal and relativistic particles resonantly inter¬ 
act with magnetosonic waves via TTD on a timescale 
r -1 that is approaching the cascading timescale given by 
T kk ~ k 2 /Dkk. Here the wave-wave diffusion coefficient 
of magnetosonic modes is given by 


Dkk ~ Eph k 


(W(k)\ 


(5) 


\P V l) ' 

Thus, the cascade is suppressed for wave numbers above: 
, ^ 8i Iq ( (£i ri(M)) N " 2 


14p(E ph ) V 


( 6 ) 


D - M = Br P (^)' , L wmdt ’ (3) 

where /3 is the thermal-to-magnetic pressure ratio, and c 
is the speed of light. The energy density IF of a mode 
in a magnetized plasma stems from both electromagnetic 
fields and resonant particles. For a high-/3 plasma, the 
pitch angle averaged ratio of beta-weighted magnetic-to- 
total energy density saturates to ((3\Bk\ 2 /2W) « 10 1 ’ 4 
(see figure 2 in Q). The pitch angle of the CR momentum 
with the magnetic field orientation is given by 6 , and 


where 2n/k cu t ~ 0.1 — 1 kpc in the ICM. This con¬ 
stitutes an effective mean free path for CRs, unless 
plasma instabilities can mediate interactions between 
turbulence and particles on smaller scales 12]. In 
this work we only consider damping via TTD due to 
thermal electrons, and neglect subdominant damping 
with thermal protons and relativistic particles. The 
latter will be subdominant in the ICM for a CR to 
thermal energy density ratio < 10% Q, which is al¬ 
ways satisfied. The azimuthally averaged turbulent 
damping rate from thermal electrons [ 8 [ in a high-/3 
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plasma is (r e ) ~ (kV ] ph \/ 37 ra;/ 20 cxp(—5a;/3) sin 2 6) rs 
0.0435fcV' p h, where x = (m e /m p )/cos 2 6. The magnetic 
field in the ICM is typically ~ pG. To compute the syn¬ 
chrotron surface brightness profiles, we use the profile of 
the magnetic field strength derived from Faraday rota¬ 
tion observations of Coma [32j in combination with the 
density profile derived from X-ray measurements 


Results and Discussion. Let us first consider the 
two models which rely on secondary electrons. After 
turbulent reacceleration, the volume-weighted, relative 
CRp energy density and CRp number density inside the 
RH for M-turbulence (M-streaming), are found to be 3 
( 2 ) % and 3.0 x 10 ~ 8 (4.5 x 10“ 8 ), respectively. As we 
will see later, these densities are just of the right or¬ 
der of magnitude to reproduce radio observations in the 
Coma cluster. In addition we predict the gamma-ray flux 
within the virial radius of the Coma cluster from CRps 
that produce decaying neutral pions for M-turbulence ( M- 
streaming) with F' 7 (> 500MeV) = 1.6 x 10 _lo (2.3 x 
10 _1 °) phs _ 1 cm -2 . Both fluxes are well below current 
limits set by Fernh-LAT (2(|, and will be challenging to 
probe in the near future. The spectral index of the CRp 
distribution is relatively steep a p ~ 2.6 for the CRp en¬ 
ergies E > 10 GeV that are relevant for the injection of 
radio-emitting secondary CRes. The steep spectrum is 
ultimately a consequence of our test particle model for 
Fermi-I acceleration [ij|, where we steepen the spectral 
index to avoid acceleration efficiencies above £ p = 15%. 

In Fig. [l] we find that all three scenarios in which the 
seeds undergo Fermi-II reacceleration can reproduce the 
Coma RH profile at 352 MHz. In the panel Brunetti et 
al. (2012) we show that without CR streaming or a flat 
turbulent profile, our simulations of reaccelerated CRs 
produce radio profiles that are too steep. Indeed, even 
using the assumptions of previous work - where complete 
freedom in the seed population was allowed it is not 
possible to reproduce observations in both frequencies. 
This signals that this problem is generic and requires ei¬ 
ther additional modifications to the plasma physics of 
Fermi-II acceleration or a better understanding of poten¬ 
tial observational systematics. 

In principle, reacceleration via TTD leads to spectral 
steepening with particle energy due to the inefficiency of 
the acceleration process to counter the stronger cooling 
losses with increasing energy. Since synchrotron emission 
peaks at frequency v s yn ~ 1 H//xG( 7 / 10 4 ) 2 GHz, this 
translates into a spectral steepening of the total radio 
spectrum (see Fig. [2|). A given radio window samples 
higher energy electrons for a decreasing field strength 
in the cluster outskirts. Hence, the spectral steepen¬ 
ing with en erg y should translate into a radial spectral 
steepening [13]. However, because of the weak depen¬ 
dence of the electron Lorentz factor on emission fre¬ 
quency (7 oc y/Vsyn), this effect is only visible in our 
simulations for u syn > 5 GHz. Most importantly, our 
simulated fluid elements at a given radius sample a broad 


distribution of shock history, density and temperature, 
which implies very similar synchrotron brightness pro¬ 
files at v syn = 352 MHz and 1.4 GHz. The discrepancy 
of the observed and simulated 1.4 GHz profiles could in¬ 
stead be due to systematic flux calibration error in single 
dish observations. Interestingly, we can match the 1.4 
GHz data if we reduce the zero point by adding 10% of 
the central flux to every data point; this flattens the outer 
profile [13] • Alternatively, this may point to weaknesses 
in the theoretical modeling of the particle acceleration 
process and may require a stronger cutoff in the particle 
energy spectrum. 

In Fig. [5] we show that our three models that include 
Fermi-II reacceleration can individually reproduce the 
convexly curved total radio spectrum found in the Coma 
cluster. Seed CRs that do not experience turbulent reac¬ 
celeration have a power-law spectrum in disagreement 
with observations. In order to match both the spatial 
and spectral profiles in Coma, we can constrain the ac¬ 
celeration efficiency for the strongest shocks in our three 
models M-streaming , M-turbulence , and M-primaries to 
C P < 0.15, £p < 0.05, and Ce < 0.004, respectively. Fol¬ 
lowing the Mach number (AJ)-dependence of the accel¬ 
eration efficiency suggested in Q, the efficiency in weak 
shocks (At ~ 2.5— 3.5) that dominates the CR distribu¬ 
tion function, has an acceleration efficiency for protons 
(j p ~ 0.0001 — 0 . 01 , and for electrons ( e ~ 0 . 001 . 

Conclusions. The standard reacceleration model 
for RHs requires a population of seed electrons to un¬ 
dergo turbulent reacceleration. These seeds are generally 
thought to be secondary electrons from hadronic CRp 
interactions. In this work we use cosmological simula¬ 
tions to derive a population of seed CRps originating 
from structure formation shocks and merger shocks dur¬ 
ing the cluster build up. The resulting secondary popu¬ 
lation is inconsistent with RH observations. We propose 
3 possible solutions that all produce gamma-ray emission 
below current upper limits and that reproduce both the 
spectrum and the surface brightness profiles of the Coma 
RH: (i) injected turbulence that is flatter than in pre¬ 
viously adopted models, (ii) streaming CRs, (iii) shock 
accelerated CR electrons injected with K ep ~ 0.1. We 
will pursue further implications and distinguishing char¬ 
acteristics of these competing models in future work. 

Acknowledgments. We thank Josh Wiener for dis¬ 
cussions on CR streaming. We are also grateful to 
Lawrence Rudnick for discussion on uncertainties in the 
1.4 GHz radio data. A.P. is grateful to the Swedish re¬ 
search council for financial support. S.P.O. thanks NASA 
grant NNX12AG73G for support. C.P. gratefully ac¬ 
knowledges financial support of the Klaus Tschira Foun¬ 
dation. 


* Electronic address: apinzke@fysik.su.se (AP) 
t Electronic address: peng@physics.ucsb.edu (PO) 







5 


M-turbulence: 


M-streaming: 



M-primaries: 



Brunetti et al. (2012): 




FIG. 1: Radio surface brightness profiles of Fermi-II reaccelerated CR electrons of a simulated post-merging cluster similar to 
Coma. We compare profiles at 352 MHz (blue lines and crosses |34j]) to those at 1.4 GHz (green lines and crosses [HI). The red 
crosses show the reprocessed 1.4 GHz data, where a zero level of about 10 % of the central value is adopted. The solid lines show 
predicted emission from a reaccelerated fossil population, while dotted lines show emission from a fossil population without 
reacceleration. The panels show the emission from CR protons and secondary electrons reaccelerated by a flat turbulent profile 
(upper left panel), secondary electrons generated by streamed CR protons (upper right panel), primary electrons (lower left 
panel), and simulated secondary electrons together with previous estimates [13| for the Coma cluster (lower right panel). 


* Electronic address: christoph.pfrommer@h-its.org (CP) 

[1] G. Brunetti and T. W. Jones, International Journal of 
Modern Physics D 23, 1430007 (2014). 

[2] F. Govoni, T. A. Enfilin, L. Feretti, and G. Giovannini, 
A&A 369, 441 (2001). 

[3] L. Feretti, G. Giovannini, F. Govoni, and M. Murgia, 
AAPR 20, 54 (2012). 

[4] Planck Collaboration, P. A. R. Ade, N. Aghanim, M. Ar- 
naud, M. Ashdown, F. Atrio-Barandela, J. Aumont, 
C. Baccigalupi, A. Balbi, A. J. Banday, et al., A&A 554, 
A140 (2013). 

[5] F. Miniati, D. Ryu, H. Kang, and T. W. Jones, Astro- 
phys. J. 559, 59 (2001). 

[6] C. Pfrommer, MNRAS 385, 1242 (2008). 

[7] G. Brunetti, G. Setti, L. Feretti, and G. Giovannini, MN¬ 
RAS 320, 365 (2001). 


[8] G. Brunetti and A. Lazarian, MNRAS 378, 245 (2007). 

[9] V. Petrosian and W. E. East, Astrophys. J. 682, 175 
(2008). 

[10] D. O. Chernyshov, V. A. Dogiel, and C. M. Ko, Astro¬ 
phys. J. 759, 113 (2012). 

[11] C. L. Sarazin, Astrophys. J. 520, 529 (1999). 

[12] G. Brunetti and A. Lazarian, MNRAS 412, 817 (2011). 

[13] G. Brunetti, P. Blasi, O. Reimer, L. Rudnick, 
A. Bonafede, and S. Brown, MNRAS 426, 956 (2012). 

[14] A. Pinzke and C. Pfrommer, MNRAS 409, 449 (2010). 

[15] F. Vazza, C. Gheller, and M. Briiggen, MNRAS 439, 
2662 (2014). 

[16] T. Enfilin, C. Pfrommer, F. Miniati, and K. Subrama- 
nian, A&A 527, A99+ (2011). 

[17] J. Wiener, S. P. Oh, and F. Guo, MNRAS 434, 2209 
(2013). 


























6 



V [Hz] 


FIG. 2: Radio synchrotron flux as a function of frequency. 
Red lines are derived from simulations, while the black crosses 
are compiled from observations [36]. The lines show the emis¬ 
sion from secondary electrons generated by streamed CR pro¬ 
tons (solid line), CR protons and secondary electrons reaccel¬ 
erated by a flat turbulent profile (dashed line), and primary 
electrons (dash-dotted line). We contrast the reaccelerated 
populations to a population of DSA accelerated CR protons 
that produce secondaries (dotted line). 
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